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Abstract — Orthogonal designs are fundamental mathematical 
notions used in the construction of space time block codes for 
wireless transmissions. Designs have two important parameters, 
the rate and the decoding delay; the main problem of the theory 
is to construct designs maximizing the rate and minimizing the 
decoding delay. 

All known constructions of CODs are inductive or algorithmic. 
In this paper, we present an explicit construction of optimal 
CODs. We do not apply recurrent procedures and do calculate 
the matrix elements directly. Our formula is based on a cubic 
function in two binary n-vectors. In our previous work (Comm. 
Math. Phys., 2010, and J. Pure and Appl. Algebra, 2011), we 
used this function to define a series of non-associative algebras 
generalizing the classical algebra of octonions and to obtain sum 
of squares identities of Hurwitz-Radon type. 

Index Terms — Orthogonal designs, decoding delay, maximal 
rate, peak-to-average power ratio, space-time codes, generalized 
octonions. 



I. Introduction 

Orthogonal designs first appeared in the classical work 
of Hurwitz |13|, |14| and Radon fT7|, in order to solve 
the problem of sum of squares identities (also known as 
composition of quadratic forms). This problem can be for- 
mulated in different ways and related to many mathematical 
questions (normed division algebras, vector fields on spheres, 
Clifford modules, immersion of projective spaces in euclidean 
spaces...) arising in different fields. The general problem is 
widely open and keep inspiring work of many mathematicians, 
see ifTSl and ll20l for surveys. In the 1970's, orthogonal designs 
and their generalizations have been extensively studied from 
combinatorial and number theoretic viewpoints, see Geramita 
et al. |6|-[9| and references therein. 

Orthogonal designs keep attracting much attention, since 
they are used to construct space-time block codes for wireless 
communication with multiple transmit antennas. This idea was 
introduced by Tarokh, Jafarkhani and Calderbank [.12] . as a 
generalization of the Alamouti scheme (Sj for wireless com- 
munication with two antennas. Space-time block codes built 
out of the orthogonal designs achieve full transmit diversity 
and have a simple maximum likelihood decoding algorithm. 

In this paper, we describe a method of construction of 
orthogonal designs. Unlike all known constructions which are 
inductive, i.e., use block matrices of small sizes to construct 
bigger matrices, our construction calculates elements of the 
matrices directly. In particular, we construct designs satisfying 
optimal criteria of ifTOl and ||T], ||2l. We also construct designs 
of type fTSl and [l5| defined by matrices that have no zero 
elements. 
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A. Definitions and known results 

Definition 1: A real orthogonal design (ROD) of type 
[p, n, k] is a matrix G of size p x n with real entries 

0, ±xi, ■ ■ ■ , ±Xk, satisfying 



G^G= (xj 



^n 5 



where G^ is the transpose matrix of G. 

Definition 2: A complex orthogonal design (COD) of pa- 
rameters [p, 71, k] is a matrix H of size p x n with complex 



entries 0, ±zi, 
satisfying 



, izfe, and their conjugates ±z*, 



,±zt 



where H* is the complex conjugate transpose of H. 

Definition 3: Given a [p,n,k]-(R or C)OD, the ratio - is 
called the rate of the design and the parameter p is called the 
decoding delay of the design. 

The main problem in the construction of real or complex 
orthogonal designs is to maximize the rate - and minimize 
the delay p for a given n. The following answers have been 
provided 

1) [p,n, fc]-ROD of rate 1 exist for all n, and in this case 
the minimum delay is p — 2*'^"^ where 



Sin) 



ra-l 
2 

K+1 

2 



V 2 ^ 



ifn = 2,4,6 mods 
if n = 1, 7 mod 8 , 
if n = 3, 5 mod 8 , 
if 71 = mod 8. 



This is a way to formulate the classical theorem or 
Hurwitz and Radon. 

2) Using a doubling process of ROD of rate 1, Tarokh et al. 
L12| obtain COD of rate i and decoding delay 2*(")+i; 
We will denote by TJC„ this class of CODs, the 
parameters are 

2'5(")+i jj 2''^"^ 

3) Liang fTOl proves that the maximal rate of a [p, n, k]- 
COD with 71 7^ p is i + i, if 71 is even, and i + ;i+T' 
if 71 is odd. 

4) Adams et al. H] and 121 find a tight lower bound for 
the decoding delay p in a non-square COD achieving 
the maximal rate given by a binomial coefficient. Let 
71 = 2771 — 1 or 77 = 2771, then 



P 



> 



2m 

771—1 



for 71 = 0, 1, 3 mod 4 and 



p>2 



2m 

?71 — 1/ ' 



for n = 2 mod 4. We will denote by LA„ the class of 
CODs achieving the maximal rate and minimal decoding 
delay. 

5) Liang [10] and Lu et al. ITTI give algorithms to produce 
designs of type LA„. Another construction is given 
in |4|. 

6) Das and Rajan Q construct CODs of rate i and 
decoding delay 2*^"'' . We will denote by DR„ this class 
of CODs, the parameters are 



2'5(n) n, 2'5(")-i 



It is interesting that the decoding delay of these CODs 
is twice lower than that of TJC„. 

Let us stress that CODs of rate i are of interest, since for 
large values of n the maximal rate is almost ^ For instance, 
for n — 12 the COD LA„ has rate j^ ^^'^ decoding delay 
792, whereas DR„ has rate ^ and decoding delay 64, see (5] 
for a comparative table. 



B. Main results and organization of the paper 

The goal of this article is to present a unified construction 
of RODs and CODs. 

We start with an explicit formulas for the matrices G of 
size 2"^ X 2*" satisfying the conditions of ROD. The rows and 
columns of the matrices are labelled by elements in the set Zj, 
i.e., of r-vectors with coefficients and 1. The entries in the 
matrices are given by an explicit function / : Zg x Zj — J> Z2. 
We then explain an easy way to reduce a RODs to a CODs. 
With our method we construct RODs with parameters; 
• [2*(2r)^2r, 2*(2'-)], 

that will produce CODs with same parameters as TJC„ 
(using a doubling process) and with same parameters 
as DR„ (using a reduction process), provided n ^ 1 
mod 8; 



2((:,t\),2r,2C) 



4(,'-+\),2r, 2(M 



, if r = 0, 1, 3 mod 4, 



if r = mod 4, 



where m is defined by r = 2m — 1 or 2m. These RODs 
will reduce to COD with same parameters as LA„. 

The paper is organized as follows. The next section contains 
the main ingredients of our approach. We construct RODs with 
parameters that are twice the parameters of the optimal CODs. 

Section Hn] describes the reduction from the [p, n, fc]-RODs 
to [f , |, |]-CODs, leading to optimal CODs of type LA„. 

Section |IV] presents a procedure that allows us to construct 
a [f , n, fc]-COD out of a [p, n, fc]-ROD, provided the ROD is 
stable under the duality. We thus obtain the CODs of types 
TJC„ and DR„. Let us mention that the corresponding 
matrices have no zero entries. 

Proofs of technical statements, as well as properties of the 
binary functions we use, are collected in Appendix. 



11. General construction of RODs 

A. Combinatorics over Z2 

We denote by K the set of r-vectors u = {ui, . . . ,Ur), 
where u^ = or ][^ The Hamming weight \u\ of an element 
is the number of non-zero component, i.e. 

\u\ -=#{«,: = l}i<,<^- 

The sum of two elements u and v is just the sum componen- 
twise modulo 2. Every element is a sum of the basis vectors 

£^ = (0,. .. ,0,1,0,. ..,0), 

with 1 at j-th position. We will also consider the element of 
maximal weight r: 

£=(1,1,. ..,1). 

We will use the involution on Z2, that we call the "hat duality": 



u = u + e , 



(1) 



i.e., the change of 1st coordinate. 

The following function in two arguments / : Z2 x Zj — 5- Z2 
plays the key role in our approach: 



I{U,V)= ^ {UjUjVk + UjVjUk + VjUjUk) + y^ UjVj, 
i<j<k i<j 



We will also use the function in one variable a{u) := f{u, u), 
given explicitly by 

a{u) ~ 2, UiUjUk + 2, UiUj. 

i<j<k i<j 

The value of a{u) depends only on the weight of u: 

if |u| = mod 4, 
1 , otherwise. 



a{u) 



The function / is used in all the constructions to determine 
signs, while a is used as a "statistic" to select good elements 
of Z2. Properties of / and a are presented in Appendix. 

B. General construction of RODs 

In this section, we construct (p x n)-matrices whose rows 
and columns are indexed by subsets W dH^ and F C Zj of 
cardinality p and n, respectively. 

We define the matrix Gu, m e Zj by 



Gu — 



■■■ Gl 



\ ■■ J 

where the entry in position (w, v) is 

'' (-l)/K^), if U^V + W, 

0, otherwise. 






'We use the notation Z2 = {0, 1} for the abeUan group of rank 2, other 
notations: F2 and Z/2Z are also often used. 



The following properties are obvious: 

1) the matrix Gu has at most one nonzero element on each 
row and on each column; 

2) ifV^W = Zl, then G„ is a (2'' x 2'')-matrix with 
exactly one non-zero element on each row and column. 

Definition 4: We call [/, V, W an admissible triple if the 
following two conditions are satisfied: 

1) 

w ^U + v, 

i.e., u + V (i W for all w G [/ and v ^ V and every 
element w (iW can be written in the form w = u + v. 
2) If a non-zero element w ^W decomposes in two ways: 
w — u + v^u' + v' then 



a{u + u') = a{v + v') = 1. 



(2) 



Theorem 1: IfU,V^W is an admissible triple, then one has: 

(i) Gl Gu = In, for all u€U. 

(ii) Gl Gu' + Gl, Gu = 0, for all u^u' eU. 

This theorem is proved in |16| and |,15j . For the sake of 
completeness, we include the proof into Appendix. 

Corollary 1: If U, V, W is an admissible triple, then 
(i) the matrix 

G — y ^ XuGu 

ueu 

is a ROD with parameters [#W, #V, #f7] in real variables a;„, 
u ^ U, where # is the cardinality of a set; 

(ii) in the case of rate 1, i.e., where k — p, the matrix G 
has no zero entries. 

Our next task is to construct admissible triples [/, V, W . 
We will use the fact that the function a vanishes only on the 
elements whose weight is a multiple of 4. Note that the easiest 
way to guarantee condition (|2]l is to choose the set V so that 

a{v + v') = l for all w, v' G V. 

C. RODs of rate 1 

In this section, we provide triples of sets U, V, W that 
produce real orthogonal designs 



G — y ^ XuGu 



ueu 

of rate 1, with minimum delay, i.e. the parameters of G are 

[2^'^"\n, 2''(")] , (provided n ^ 1 mod 8). 

Case r = 0, 1, 2 mod 4. One chooses the following subsets 

V = {£^e'^, l<i<r}, 

U = W = 1.1, 
where "is the duaUty dTJ. Then G is a [2'^, 2r, 2'']-ROD.E| 



- This ROD is optimal, except for the case r = mod 4, where, according 
to the Hurwitz-Radon theorem, there is a [2"", 2r + 1, 2'']-ROD. We do not 
dwell here on a more involved construction to produce such a ROD. 



Case r = 3 mod 4. One chooses the following subsets 

V = je, £, e-J', P, 1 < j < r|, 
U = W = 11, 
the G is a p'', 2r + 2, 2''] -ROD. 

D. Non-square RODs of rate i + — !- 

All the RODs below have maximal rate i + tt-, when 
r — 2m or 2m — 1. 

Case r = 1, 2 mod 4. Consider r = 2m, — 1 or r = 2m and 
choose the set U of the elements of weight m and their dual, 
the set V is chosen as in the first case: 

U — {u, u : \u\ — m} , 

V = [e^, d, l<i<r}, 

It follows that the space W = U + V \^: 

W = {w : |u| = rn — 1, 777,, 777, + 1} y 
{u : Ui — 1 , |u| = 777 + 2} y 

{m : Ml = , |7i| = 77T, — 2} . 

The matrix G is a ROD with parameters 



'r+l 



i)>2r,2(;j 



for odd r and even r, respectively. 



-i)>2r,2(,;) 



Case r — {) mod 4. Consider r — 2m (where 771, is even) and 
choose the following subsets 



U 

V 

w 

then G is a 



{u : Ui = 1 , \u\ = 777,} y 

{u : til = , juj = 777, — 1} , 



}• 



= <e,e, e^, e^, 2 < j <n 

= {u : Ul = 1 , juj = 777 — 1, 777 + 

{u : ui — , |?i| = 777 — 2, 777} 



i}U 



2( \),2r,2('-\) 



-ROD. 



Case r — 3 mod 4. Consider r' :— r + l and apply 
the previous case with r' = mod 4 to obtain a 
-ROD, where r = 2m — 1. Removing 



2r+\),2r + 2,2('" 



two columns, we obtain a 



2r+\),2r,2r) 

\m—l/ ' ' \mJ 



-ROD. 



In each of the above cases, condition (|2|i is satisfied for all 

v,v' e V. 

To finish this section, let us mention that the binary nu- 
meration have already been efficiently used in l|4], ||5] to 
construct RODs and CODs of maximal rate. In particular, 
subsets of Z2 similar to our sets U, V and W were described. 
The main difference of our approach is the function / and 
explicit construction of the matrices. 



III. Reduction from ROD to COD 

In this section, we present a procedure to reduce a 
[2p, 2n, 2fc]-ROD to a [p, n, fc]-COD. Such a procedure is not 
always possible, it requires nice properties of sets U, V, W. 
We first describe the general procedure of reduction and then 
apply it to RODs of rate 1 constructed in Section III-CI 

A. The general procedure 

The main idea is to use a duality 

^ : Z'2 ^> IT2 



and 



Z2 satisfies 



defined by u = u + e, where the element e G 

/(e, e) = 1, and to choose sets U, V and W stable under the 

duality: 

U ^U, V = v, w = w, 

In practice, we use the hat duality ([T]i, i-C-, e = e^. 

Given a ROD of type [2p,2n,2k] defined by sets U,V 
and W in Z2, our goal is to reduce it to a COD with 
parameters [p,n,k]. The method consists in two steps. 
First, we introduce splitting of the sets U, V, W, in order to 
decompose the matrices G = X^uec/ ^uGu into admissible 
(2 X 2)-blocks. Second, we replace the admissible blocks by 
complex variables, 2„ ~ Xu + ixu or z* = a;„ — ixu. 



STEP 1: We fix the following sphtting of U 

Uo :- {ueC/:/(u,e) = 0}, 
Ui := {ueU:fiu,e) = l}. 



(3) 



Note that Uq = Ui, cf. Property (b) of / in Appendix. 

We now need to find subsets Vq, Vi, Wq and Wi satisfying 
the following conditions 



(4) 



where |J denotes the disjoint union. 

These splitting induce a natural decomposition of the ma- 
trices Gu into (2 X 2) -blocks whose columns are labelled by 

{v, v) eVo X Vi and the rows by {w, w) G Wq x Wi : 



V = 


-- ^oUl^i, 


Vi = ^o; 


w = 


= WoUWi, 


14^1 = 14^0 


^0 = 


-- Uo + Vo 


= Ui + Vu 


A^i = 


-- Uo + Vi 


= Ui + Vo, 



( 



V V 



Gu = 



\ 



G-u 
1, 



W 
W 



V / 

where u = v + w and sou = v + w — v + w. 

For u ^ Uq (and therefore u G C/i), the non-zero blocks are 
of the form 

■ (_i)/K^') \ 

(_!)/(",«) J 



/^W,V 



g:. 



(_i)/('i,«) 




non-zero blocks are located at the same place in Gu and Gu. 
Moreover, since / is linear in the 2nd variable, 

f{u,v) ^ f{u,v) + f{u,e) == f{u,v), 

f{u,v) = f{u,v) + f{u,e) = f{u,v) + l, 

so that the entries in the blocks of G„ are of the same sign 
and those of Gu are of the opposite sign. 

STEP 2: The matrix G = J2ueu ^uGu decomposes into 
(2 X 2)-blocks, and the non-zero blocks are of two types 



(Tl) ± 









or (T2) ± 



Xu 



We construct a complex matrix H from G by substituting to 
the block (Tl) the complex variable ±2;„ and to the block (T2) 
the complex conjugate variable ±z*. More precisely, the entry 
of H in position (w, v) G Wq x Vq is 

^„,, ^ f (-l)/(''+--'')z„+^, if /(^M^, v) = f{v + w, v), 

\ (-l)^^''+"'''')<+^, otherwise, 
if u + w G Uq, and H^'"" = 0, otherwise. 

Theorem 2: The constructed matrix H defines a COD with 
parameters [p, n, k]. 

B. CODs of parameters LA„ 

As application of the above procedure, let us reduce the 
RODs constructed in Section III-DI in order to obtain the 
optimal CODs of type LA„. We need to describe here the 
subsets Uq,Ui,Vq,Vi,Wq, Wi satisfying Q and dUi. 

From the expression of / we see that /(u, e^) depends only 
on the class \u\ mod 4. More precisely 



\u\ mod 4 





1 


2 


3 


f{u,e') if wi=0 








1 


1 


/(w,e^) if ui = l 





1 


1 






for an arbitrary w G ^3. 

Case r = 1, 2 mod 4. Let now m G C/, so that |u| — m and 
r = 2m — 1 or 2m. In this case, m, is necessarily odd. 
• if TO = 1 mod 4, then for u ^ U we have 

f{u,e^) = 0^=^ui=0, 
in other words, 

C/o = {ugC/|ui =0}, Ui = {ueU\ui = l}. 
We easily check that 
Vo = {veV\vi = 0}, Vi^{veV\vi^l}, 
Wo^{w eW\wi=0}, Wi^{w(zW\wi = l}. 
satisfies property (|4}. 



• if ?7i = 3 mod 4, then for u E U we have 

in other words, 

Uo = {ueU\ui = l}, Ui^{ueU\ui^Q}. 

We easily check that 
Vo = {veV\vi = 0}, Vi^{veV\vi^l}, 
Wo^{w eW\wi^l} Wi^{weW\wi=0}. 

again satisfies property (|4]l. 

Case r = mod 4. In this case, m is defined by r = 2m, so 
that m is even. 

• If TO = mod 4, then for u E U we have 

in other words, 

UQ^{ueU\ui^l}, Ui^{u£U\ui^O}. 
We easily check that 
Vo^{veV\vi^O}, Vi^{veV\vi = l}, 
Wo = {w eW\wi^l}, Wi:^{w£W\wi^O}. 

have the desired property. 

• If TO, = 2 mod 4, then for u £U we have 

/(u,ei) = 0^=^ui = 0, 
in other words, 

[/o - {u e C/ 1 ui - 0} , C/i = {u e C/ I ui = 1} . 
We easily check that 
Vq = {v£V\vi = Q}, Vi = {veV\vi = l}, 
Wg^{wC,W\wi^Q} Wi = {weW\wi^l}. 
have the desired property. 

IV. Induction 

The second procedure that we call induction allows one to 
transform a [p, n, fc]-ROD to a [p,n, |]-COD. This induction 
consists in complexification composed with reduction; it can 
be applied whenever the set U, V and W are stable under an 
involution u = m + e for an element e satisfying /(e, e) = 1. 

A. The general construction 

STEP 1 : We consider the splitting of [/ as in (O and define 
the following subsets of Z^+^ = Z^ x U Z^ x 1 

U' = C/oxOUL/ixl, 

F' = T/xOUFxl, 

W = W^xOUM^xl. 

For each u E U, we embed the previous matrix G„ into a 
twice bigger matrix, G(„ ,-), where (u, t) G t/', with columns 



indexed by W and rows indexed by V. The matrix Gi^^.t) is 
composed by (2 x 2) -blocks: 



( \ 



Cr(u,T) 






\ 






/ 



where the non zero blocks correspond to {u, r) = (u + w, r) 
and coincide with those of Gu- These data provide a 

[2p,2n,k]-ROD. 

STEP 2: For each u e Uq, we define 

G^ ~ XuG(u,o) + 2;iiG(„ 1). 

These matrices decompose into blocks that are all of type 
(Tl) or (T2). We then apply the reduction procedure to obtain 
a [j5,n, |]-ROD. 

B. CODs with parameters TJC„ and DR„ 

Consider the RODs of rate 1 constructed in Section III-CI 



Following 1 12 1, one can apply the following obvious dou- 
bling process: 

(2) ._ ( G 

a 



^(2 



where the variables x^ in G are now considered as complex 
variables, and G' is copy of G associated to the conjugate 
variables, i.e. is defined by 



G — 2_^ x'^Gy 



ueu 
The matrix G^^-* is a COD of rate i, with parameters 

that are precisely the parameters of TJC„, for even n. 
Removing a column, we are led to CODs with parameters 
TJC„, for odd n, except for n = 1 mod 8. 

Applying the induction procedure to the RODs of rate 1, 
leads to CODs with parameters 



2r)-l 



25(2r)^ 2r, 2^( 



that are precisely the parameters of DR„, for even n. Again, 
removing a column, we obtain CODs of type DR„, for 
odd n, except for n = 1 mod 8. 

In both cases, the obtained CODs have no zero entries. 

The missing case n — I mod 8 escapes from the technique 
used in this paper This is due to the fact that we do not obtain 
a [2'', 2r + 1, 2''] -ROD with n = mod 4. 



V. Appendices 
A. Properties of the functions f and a. 

The function / has quite remarkable properties that we 
briefly discuss here. 

It is impossible to reconstruct the function / from the 
function in one variable a (which is nothing but the restriction 
of / to the diagonal in Zj x Zj). However, a contains the 
essential characteristics of /, such as its symmetrization. 

1) First polarization formula: 

f(u, v) + f{v, u) = a{u + v) + a{u) + a[v). 

2) Second polarization formula: 

/(u, v) + f{u, V + w) + f{u + V, w) + f{v, w) = 

a{u + V + w) 

+a(u + v) + a{u + w) + a{v + w) 

+a(u) + a(v) + a{w). 

These properties can be checked directly. Note that the expres- 
sion in the right-hand-side of [Til is called the coboundary of a, 
it has a deep cohomological meaning. The expression in the 
left-hand-side of |2]l is the coboundary of /, its measures the 
non-associativity of a certain algebra defined by /, see IfT^ . 
Finally, the expression in the right-hand-side of |2]i is called the 
polarization of the cubic form a. Let us mention that, unlike 
the theory of quadratic forms, the theory of cubic forms is not 
well developed in characteristic 2, not much is known. 

Let us also give here more elementary properties of / 
already used in the above constructions: 

(a) Linearity of / in 2nd variable: 

f{u,v + v')^f{u,v) + !{u,v'). 

(b) Pseudo-Unearity in 1st variable: 

f{u -^v,v) = f{u, v) + f{v, v). 



We invite the reader to consult lil6J for more information 

about / and a. 

B. Proof of Theorem [7] 

We apply the formula for matrices multiplication. The 
coefficient in position {v,v') in the product G^Gu' is 



{GlGu'T^ 



(_l)/(«,")+/(«'y) ifv + v' = u + u', 







otherwise. 



Proof: Rewrite the left-hand-side using v' = u+u'+v and 
the linearity in the 2nd variable, after cancellation of double 
terms one obtains 

f{u', u) + f{u\ u') + f{u, u) + f{u, u'). 

This reduces to a{u + u') using the first polarization formula. 

■ 
Theorem [T] follows. 

C. Proof of Theorem |2] 

First notice that in each column of the matrix H the symbol 
Zu appears exactly once ("symbol z„" means one of the 
following four elements: ±z„,±z*). This implies that the 
diagonal entries in H*H are all equal to 



u&Uo 



w? 



It remains to show that the non-diagonal entries in H*H are 
all zero. We show that, in the hermitian product of two distinct 
columns of H, the terms pairwise cancel. 

Consider the four entries of H, in position (w,w), {w' ^v), 

{w',v) and {w',v'). 



( 



H 



\ 



z\ 



Z3 



Z2 



Zi 



V 



w 



w 



I 



Case I: there exist m, u' in t/o such that 

M = W + W = ti'+ w', u' = W + w' — w' + W. 

In this case, the four entries are non zero and one has 

Zl,Z4 e {±Z„, ±Z*^ , Z2,Zz e {±Zu', ±2*'}- 

The corresponding blocks in the matrix G 

( : : \ 

... Ai • 



This impHes that G^G^, = /„, and G^G^, + G^^.G^ = if 
and only if 

(■_!')/(". '") + /("'^"') _|_ ('_n/('"':'") + /(«:'"') — 

whenever t; + w' = u + it! . The above condition is equivalent 
to 

/(m, v) + /K, v') + /K, v) + /(m, v') = 1. 

Lemma 1: \f u + u' — v + v' then 
f{u, v) + fiu', v') + fiu', v) + fiu, v') = a{u + u'). 



G = 



come from x^Gu 
satisfy 

This translates to 



A. 



Ao 



A, 



I 



XuGu + Xu'Gu' + x-.G;-. and therefore 



AlA2 



AlAi = 0. 



ZJJ'Z2 + 23^4 = 0. 

Case II: there do not exist u, u' in C/q such that 

u^v + w^v' + w' , u' = V + w' — v' + w. 
In this case at least one of the following situations holds 

zi = Zi — or 22 = Z3 = 0, 



and, again, zj'z2 + z|z4 = 0. 

We have proved that the columns of H are pairwise orthog- 
onal (with respect to the hermitian product). And we conclude 
finally that 

H*H = [Y,u£Uo\^u\ ) In 

where n = #Vb. 

D. Orthogonal designs and Hurwitz problem of sums of 
squares 

It is well-known that the existence of [p, n, k]-ROD is 
related to the Hurwitz problem on composition of quadratic 
forms, IIT4I . IIT7I (see also 1201 for a survey). 

Definition 5: A Hurwitz sum of squares identity (SSI) of 
size [p, n, k] is an identity 



(«? 



(bl 



bl) 



(5) 



where ct are bilinear expressions in a^ and bi with integral 
coefficients (the elements ai,bi,Ci's are considered here as 
real variables). Such an identity will be referred as a [p, n, k]- 
identity. 

It is known that if such an identity holds then the integral 
coefficients in the expressions of q's can be chosen among 
{0,1,-1}. Hurwitz proved the following fundamental theo- 
rem. There exists a [p, n, k]-ROD if and only if there exists a 
[p, n, k]-SSI. 

Let us recall here how the equivalence can be established. 
Since c's are linear in a's and 6's, one has 




a,Ai 



(6) 



where 6 is a column-vector with components hi and c is a 
column- vector with components Ci, and where Ai are p x n 
matrices (with entries 0,1,-1) One then easily checks that 
the identity ^ holds if and only if 



Ai A,i 



A] A, = 0, 



V J 7^ j. 



(7) 



Ei<t<fcfi*^0 is a [p, n, A:]-ROD in 



Ai Ai — J„ 



Then, the matrix A 
the variables a/s. 

A classical result of Hurwitz lfT3l states that [n, ?i,ri]-SSI 
exist if and only if n — 1,2,4,8. This statement relies on 
classification of normed division algebras, see 1191 for a survey 
on relations between division algebras and wireless communi- 
cations. The case of [n, n, fc]-SSI was solved independently by 
Hurwitz llT4i and Radon |17| , this is the origin of the famous 
Hurwitz-Radon function. 

Let us mention that our previous results |T6l and 1 031 
were formulated in terms of partial solutions to the Hurwitz 
problem. 
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